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We discuss the renormalizability of ^-derivable approximations in scalar (p 4 theory in four dimen- 
sions. The formalism leads to self-consistent equations for the 2-point and the 4-point functions 
which are plagued by ultraviolet divergences. Through a detailed analysis of the one and two-loop 
self-energy skeletons, we show that both equations can be renormalized simultaneously and deter- 
mine the corresponding counterterms. These insure the elimination of ultraviolet divergences both 
at zero and finite temperature. 
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Self-consistent, "<I>-derivable" , approximations were in- 
troduced many years ago in the context of the non- 
relativistic many body problem P, and have been 
extended to field theory (3, 0|- They have been found 
appropriate to treat systems for which the quasiparticlc 
picture is a good starting point and have recently been 
applied in this spirit to calculate equilibrium thermody- 
namics of the quark-gluon plasma [||. They are also be- 
ing used to study the dynamics of quantum fields out of 
equilibrium 

The main difficulty in implementing such approxima- 
tions in quantum field theory is their renormalization: 
from the point of view of perturbation theory, the equa- 
tions that one is led to solve effectively resum infinite sets 
of Feynman diagrams, and the existence of a procedure 
for constructing the counterterms needed to eliminate the 
corresponding divergences is not obvious. This problem 
becomes particularly acute at finite temperature: While, 
on general grounds, one expects ultraviolet divergences 
to be unaffected by the temperature (see e.g. 0), in self- 
consistent approximations temperature dependent diver- 
gences often do appear, thus casting doubts on the renor- 
malizability (see in particular Q). 

This issue has been addressed recently by van Hees 
and Knoll in a series of papers |jj 0] . The strategy put 
forward in 

.9] 

is based on an expansion of the propagator 
around the vacuum self-consistent solution, and relies on 
the real time formalism. The elimination of the diver- 
gences proceeds through the BPHZ subtraction scheme. 
This leads to a systematic and practical renormaliza- 
tion scheme where temperature dependent counterterms 
never appear. However the dissymmetrical treatment of 
the vacuum sector and the finite temperature one is un- 
satisfactory: It hides the fact that the rearrangement of 
divergences which appears to be necessary at finite tem- 
perature is also needed in most renormalization schemes 
already at zero temperature. And it does not bring out 
the specific relation between the bare and the renormal- 
ized parameters that emerges in ^-derivable approxima- 
tions. This makes it difficult, e.g., to compute the (3- 
function, or resolve the apparent discrepancy between 
the results of Refs. y| and [9j . 

We have therefore reconsidered the problem from a 
more general perspective. Our derivations use the imagi- 



nary time formalism, making the connection with conven- 
tional equilibrium field theory transparent, and allowing 
for a simultaneous treatment of the vacuum sector and 
the finite temperature one: once renormalization is done 
properly at zero temperature, the extension to finite tem- 
perature is straightforward. 

The central quantity in <I>-derivable approximations is 
$[.D], the sum of the 2-particle-irreducible "skeleton" 
diagrams, a functional of the full propagator D, which 
enters the expression of the thermodynamical potential. 
From $[.D] we may calculate the 2-point function (the 
self-energy) by functional differentiation: 



5<5>[D}/5D = ~TI. 



(1) 



This relation, together with Dyson's equation (Do de- 
notes the bare propagator): 



D- 1 = Dn 1 



U[D], 



(2) 



defines the physical propagator and self-energy in a self- 
consistent way. We shall refer to Eq. @, with II given 
by , as the "gap equation" . A further differentiation 
of $ with respect to D yields the 2-particle irreducible 
kernel 



of a Bethe-Salpeter (BS) equation 



(3) 



T(K, P) = A(K, P) 



T(K,Q)D 2 (Q)A(Q,P) (4) 



Q 



that allows the calculation of the 4-point function 
T(K, P) with a degree of accuracy comparable with that 
used in the determination of the propagator, ^-derivable 
approximations are obtained by selecting a class of skele- 
tons in 3>[D] and calculating II and T from the equations 
above. As we shall see, the renormalizability of such ap- 
proximations relies on the possibility to simultaneously 
renormalize II and T. In particular, the BS equation 
is needed to determine coupling constant counterterms 
which eliminate some divergences of the self-energy. 
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FIG. 1: The one- loop and two- loop skeleton diagrams con- 
tributing to the self-energy. These will be referred to as the 
"tadpole" and "sunset" diagrams, respectively. 

We consider in this paper a massive scalar field theory 
with a ip 4 interaction: 

£ = \ (%>o) 2 ~ ^oVo ~ jjffoVo) ( 5 ) 

and include in $ only the 2-loop and 3-loop skeletons 
(the corresponding self-energy diagrams are displayed in 
Fig. 1). This allows us to introduce the generic difficul- 
ties, deferring the systematic discussion of the general 
case to a forthcoming publication (l4|. In four dimen- 
sions, usual power counting indicates that only the 2- 
point and the 4-point functions are divergent. The di- 
vergent parts can be absorbed in local countertems cor- 
responding to a redefinition of the parameters of the 
lagrangian. We assume standard relations [3 between 
the renormalized and bare parameters: tpo — V~Zip, 
Zml = m 2 + 5m 2 , Z 2 g 2 = g 2 + Sg 2 , and SZ = Z - 1. 

The gap equation corresponding to the one-loop skele- 
ton (the "tadpole") reads: 

n = %- [ D{P) + Sm 2 (6) 
2 Jp 

where D^ 1 (P) = P 2 +m 2 +H. The notation J p stands for 
an Euclidean integral over the 4-momcntum P. At finite 
temperature, it should be understood as an integral over 
the 3-momentum together with a sum over Matsubara 
frequencies. The approximation corresponding to Eq. © 
is a simple self-consistent mean field approximation that 
has been treated many times before (see for instance 0, 
IH IH El)- We present it here in a way which will 
prepare for the more complicated two-loop example that 
we shall discuss next. 

The self-energy II is here a constant, and a single mass 
counterterm Sm 2 is in principle sufficient to eliminate the 
ultraviolet divergence. Calculating the integral in Eq. © 
in dimensional rcgularization we get: 

where p, 2 = 47re _7E /j, 2 . At this point, one could be 
tempted to absorb the whole divergence in Sm 2 , i.e., set: 

Sm 2 = ^(m 2 +U)-. (8) 

But this is not a good strategy. If, for instance, the 
calculation is done at finite temperature, II depends on 



the temperature, and so does the counterterm (JSJ), which 
we want to avoid. 

In fact, when analyzing the gap equation Eq. @ in 
terms of perturbation theory, on finds that its solution 
effectively resums an infinite set of Feynman diagrams, 
some of which contribute to the renormalization of the 
coupling constant. This is best seen by imagining solving 
this equation by iteration, a procedure which also defines 
an explicit construction of the counterterms. To do so, 
we set D(P) = D Q (P) = (P 2 + to 2 )" 1 in the r.h.s. of 
Eq. ©; one then obtains a first approximation to II on 
the l.h.s, which can then be used in the r.h.s., and so 
on. At each iteration, 5m 2 can be adjusted to absorb the 
overall divergence. But is is easy to see that, starting at 
the second iteration, a subdivergence appears correspond- 
ing to a coupling constant renormalization that needs to 
be subtracted before adjusting 5m 2 (An illustration of 
the phenomenon in the less trivial example of the sunset 
diagram is given in Fig. 2 below). New such subdiver- 
gences appear in each iteration, and to take them into 
account, a term of the form (6g 2 /2) J p D(P) should be 
added in the r.h.s. of Eq. ©. Equivalcntly g 2 should be 
replaced by g 2 = g 2 + Sg 2 in Eq. @). As we shall see, 
Sg 2 is precisely the counterterm that is needed to make 
finite the BS equation, to which we now turn. 

With, here, A = g 2 = g 2 + Sg 2 , the BS equation reads: 

r = g 2 - 9 lrj p D 2 (p), (9) 

where T is the renormalized 4-point function, and Sg 2 
is chosen so as to absorb the divergence of the integral. 
Note that this divergence does not depend on the mass 
(nor therefore on n), and for the purpose of determining 
Sg 2 we could as well use an auxiliary 4-point function Tq 
solution of Eq. © with D replaced by D . Clearly, T 
differs from T by a finite quantity only. 

We now return to the gap equation, Eq. JfjJ) with g 2 
replaced by g$, itself determined in terms of Tq by the 
BS equation, and show that its solution, II, can be made 
finite with a counterterm 5m 2 independent of II. To this 
aim, we write D = Dq + SD, where 

SD(P) = D (P)[-n]A)(P) + D r (P), (10) 

and D r (P) starts at order II 2 , so that the integral 
J p D r (P) is finite. Then, we set 

^ 2 = f J p D o( P )+Sm 2 ii, = g lj^D{P), (11) 

where only no depends on II, and when n is solution of 
the gap equation, n = Bo + B2. Next, one uses the BS 
equation to eliminate g^ in the defining equation for Bo: 

Ho = y j p 5D{P) + ^Bo J D 2 {Q). (12) 
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At this point, we have achieved our goal: while both 
integrals in Eq. (|12fl are divergent, it is easily verified 
that no divergence involves II when II is solution of the 
gap equation (so that we can replace n by II — n 2 in 
the r.h.s. of Eq. (|12f> ). The final gap equation reads: 



n 



D r (P) + U 2 1 



Q 



D 2 (Q) 



(13) 



The divergence in the last term involves only Do and can 
be absorbed in 5m 2 . For instance, in dimensional regular- 
ization with minimal subtraction Sm 2 — g 2 .m 2 /(Z2it 2 e). 
The factor multiplying II2 is To/ffg (from the BS equa- 
tion) , so that the resulting expression is indeed finite (but 
depends on the scheme). 

At finite temperature, the one-loop integral in Eq. © 
can be split into a vacuum integral and a 3-dimcnsional 
integral involving a statistical factor and giving the fol- 
lowing new contribution to II: 



n(e p ) 



(14) 



where n(e v 



l/(e 



0e v 



1) and e v = sfp 2 



n. 



Eq. (|14|) involves a temperature dependent counterterm. 
However the same manipulation as above, with IIo re- 
placed by n + fli in Eq. O, and the use of the T = 
counterterms which are calculated entirely from Do, elim- 
inate it, leaving a finite gap equation. In the mass-shell 
subtraction scheme where m is the physical mass and the 
vacuum sector is trivial (II = 0), this equation is simply: 



n(s p ) 



(15) 



Consider now the 2-loop skeleton (the "sunset"), and 
the corresponding gap equation: 

U(K) = -^J^J^ D(P)D(Q)D{K + P + Q) 

+ Ag2 ^ g2 / D(P) + Am 2 + 5m 2 + K 2 SZ. (16) 
2 Jp 

That this expression can be made finite with the in- 
dicated counterterms follows from a standard analysis: 
the counterterm Ag 2 + Sg 2 cancels the subdivergences, 
while Am 2 + Sm 2 and 5Z cancel the remaining global 
divergences (the reason behind the special writing of the 
counterterms will become clear shortly). The argument 
assumes, in agreement with Weinberg's theorem, that 
the repeated insertions of the self-energy in the propa- 
gators, as generated by iterating the gap equation, do 
not change in an essential way the asymptotic form of 
these propagators, expected to be typically of the form: 
Tl(K) ~ K 2 F(\nK) for K > m. Note that the coupling 
constant counterterm enters only a one-loop diagram: at 
this order of the skeleton expansion, there is no renor- 
malization of the vertices of the sunset diagram. Such 




FIG. 2: The sunset diagram with one sunset inserted on 
one of the propagator. The subdivergence contained in the 
dashed line box is removed by a counterterm determined by 
the Bethe-Salpeter equation. 



renormalizations would involve skeletons whose lowest 
perturbative order is g 6 . 

To proceed to the renormalization of the BS equation 
we need to take into account the fact that the self-energy 
H(K) modifies the asymptotic behavior of the propaga- 
tor, as indicated above. We then write n = II2 + IIo, 
where Tl 2 (K) is finite and contains the exact asymptotic 
behavior of H(K), and Ho(K) grows at most logarith- 
mically at large K (H 2 (K) and 5Z can be obtained by 
solving the gap equation with m = 0). We set: 



D_ 2 {K) = (K 



2 , 2 
m 



n 2 



(17) 



D_2 will play here the role of Dq in the one-loop exam- 
ple (note that D_2 takes care of field renormalization). 
Thus, we define an auxiliary 4-point function To as the 
solution to the BS equation (0} with D^ 2 as propagators. 
The equation for To(K, P) contains all the divergences of 
that for the full 4-point function T(K, P), and the renor- 
malizations of To and T involve therefore the same coun- 
terterms. To determine these, we first write the kernel of 
the BS equation as A (K,P) + Sg 2 , with: 

A (K,P) = Ag 2 -g 4 f D_ 2 {Q)D_ 2 {K + P + Q), (18) 
Jq 

and Ag 2 is chosen so as to make Ao(K,P) finite. The 
counterterm Sg 2 is then adjusted, as in Eq. ©, so as to 
eliminate the divergence of the equation: 

r (o,o) = ,5. 9 2 + Ao(0,0) 

' r T a (0,P)D 2 _ 2 (P)[Sg 2 +A (P,0% (19) 



where r (0, 0) is fixed by a renormalization condition, 
and Fo(0, P) can be obtained from the following finite 
equation : 

r (o,p) - r (o,o)=A (o,p)-A (o,o) 



r (0, Q)D 2 _ 2 (Q) {A (Q, P)-A (Q, 0)} (20) 



(A (Q,P) - A (Q,0) - 1/Q 2 for Q 2 > P 2 , so that the 
integral over Q is indeed finite.) 
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By combining Eqs. (|2*0j) and (|H)|) one gets r (0, P) in 
terms of Aq(Q,P) and Sg 2 . The result is in fact noth- 
ing but Eq. (QJ with K = 0, D replaced by £U 2 an d 
A(0, P) replaced by A o (0, P) + Sg 2 . We shall now use 
this expression of ro(0, P) to eliminate the vertex subdi- 
vergences from the gap equation. To proceed, we write 
again D = L>_ 2 + SD, where SD = £>_ 2 [-n ]-D- 2 + D r 
contains all the dependence on IIo, and call IT2 the r.h.s. 
of Eq. lfT?))l evaluated with D replaced by D- 2 . Further- 
more, we set fl (K) = (1/2) f p [A (K,P) + Sg 2 ]SD(P). 
The gap equation is II = II2 + ITo + n r where H r (K) 
is finite and goes as 1/K 2 at large K. We then write 
fi (K) = fi (K) - fi (0) + fi (K) , where fi ( K) - fi (0) 
is finite, and we express Ao(0,P) + 5g 2 in the defining 
equation for fio(0) in terms of r (0, P). We get: 



n (0) = - J T (0,P) \SD(P) 



D 2 _ 2 (P)il (P) . (21) 



For II solution of the gap equation, we can set fio = II — 
II 2 — II r in the r.h.s., and verify that the divergent terms 
linear in IIo (= II — II2) cancel, as anticipated. Using the 
resulting expression of IIo we obtain the solution of the 
gap equation for K = in the form: 



n(o) = \ J r(o, p) {n 2 (p) - fi 2 (p) - n r (P) } d 



r(o,P)D r (P) + n r (o) + fi 2 (o). 



(22) 



To isolate the remaining divergences, we write 

Sg 2 



U 2 (K)=U' 2 (K) + Sr 



Di 2 (P) (23) 



where tl' 2 (K) is finite (owing to the counterterms Am 2 
and A<? 2 ). The difference tl' 2 (K) — H 2 (K) is logarithmic 
at large K and contributes to a divergence of the first in- 
tegral of Eq. \\'2'2{i . But neither this divergence, nor those 
coming from the counterterms displayed in Eq. I|2H|I. de- 
pend on the solution no of the gap equation, and they 
can be absorbed in the mas counterterm 6m 2 . This com- 
pletes the determination of the counterterms which, as 
we have seen, can all be calculated from D- 2 . 

At this point we emphasize a special feature of $- 
derivable approximations: As we have indicated earlier, 
the renormalization of the two-loop skeleton generates a 
coupling constant counterterm for the one-loop skeleton, 
but not for its own vertices. This is a general feature, 
which persists in higher orders in the loop-expansion of <J>. 
Correspondingly, the /3-function deviates from that given 
by perturbation theory beyond the perturbative orders 
explicitly included in the skeletons [L4J. For instance, in 
the present example, the perturbative /3-function is cor- 
rectly reproduced to order g 4 (when one adds the two 
contributions of the tadpole and sunset diagrams), but 
deviates at order g 6 . 

The extension of the previous analysis to finite temper- 
ature brings no new ultraviolet difficulty. Again, we can 



separate each loop integral in the sunset into a "vacuum" 
contribution, and a contribution containing a statistical 
factor. The final expression for the self-energy takes then 
a form similar to that at zero temperature, and may be 
written as n = fi 2 + fi r + fio + fii + II3. The first con- 
tribution, fi 2 , is the same as before and does not depend 
on the temperature. The last contribution, Ha(K ), is one 
in which each of the loop integrals contains a statistical 
factor. It is finite and decreases as 1/K 2 at large K; 
thus it is not involved in any divergent term, and it can 
be regarded as a simple correction to fi r . Finally, fio is 
defined as at zero temperature, and [T^ |: 



fii(A-) 



[A (K,P) + 6g 2 ] p(p ,p) 



where p(po,p) is the spectral function of the propagator 
D, a Po denotes the sign of po, and the integral runs over 
the real po axis. As in the one-loop example, we can 
combine n! with n and show that the zero temperature 
coupling constant counterterms eliminate the apparent 
divergence depending on the temperature. 

We are grateful to H. van Hees and J. Knoll for stimu- 
lating discussions concerning their work, and to A. Reb- 
han for useful comments on this manuscript. 



[2] 

[3] 

[4] 
[5] 

[6] 

[7] 

[8] 

[9] 
[10] 

[11] 
[12] 
[13] 

[14] 



blaizot@spht.saclay.cea.fr 

iancu@spht.saclay.cea.fr 

reinosa@spht.saclay.cea.fr 

J. M. Luttinger and J. C. Ward, Phys. Rev. 118, 1417 

(1960) ; T.D. Lee and C.N. Yang, Phys. Rev. 117 22 

(1961) ; C. De Dominicis and P.C. Martin, J. Math. Phys. 
5, 14, 31 (1964). 

G. Baym and L. Kadanoff, Phys. Rev. 127 22 (1962); G. 
Baym, Phys. Rev. 127, 1391 (1962). 
J. M. Cornwall, R. Jackiw and E. Tomboulis, Phys. Rev. 
D 10, 2428 (1974). 

G. Baym and G. Grinstein, Phys. Rev. D 15, 2897 (1977). 
J. P. Blaizot, E. Iancu and A. Rebhan, Phys. Rev. Lett. 
83 (1999) 2906; Phys. Rev. D 63, 065003 (2001) 

G. Aarts, D. Ahrensmeier, R. Baier, J. Berges and J. Ser- 
reau, Phys. Rev. D 66, 045008 (2002) 

J. C. Collins, "Renormalization. An Introduction To 
Renormalization, The Renormalization Group, And The 
Operator Product Expansion," Univ. Press, Cambridge, 
UK (1984). 

E. Braaten and E. Petitgirard, Phys. Rev. D 65, 085039 
(2002) 

H. van Hees and J. Knoll, Phys. Rev. D 65, 025010 (2002) 

H. van Hees and J. Knoll, Phys. Rev. D 65, 105005 (2002) 
D. J. Gross and A. Neveu, Phys. Rev. D 10, 3235 (1974). 
L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 

I. T. Drummond, R. R. Horgan, P. V. Landshoff and 
A. Rebhan, Nucl. Phys. B 524, 579 (1998) 

J.-P. Blaizot, E. Iancu and U. Reinosa, in preparation. 



